Coherent shift of localized bound pair in Bose Hubbard model 



OO 
O 

o 

(N 

> 
O 



43 ■ 

^— > ■ 

c ■ 

=5 

cr 



> 
in 
o 

(N 



OO 

o 



X 



L. Jin, B. Chen, and Z. SongT 
Department of Physics, Nankai University, Tianjin 300071, China 

Based on the exact results obtained by Bethe ansatz, we demonstrate the existence of stable 
bound pair (BP) wave packet in Bose Hubbard model with arbitrary on-site interaction U. In 
large- U regime, it is found that an incoming single-particle (SP) can coherently pass through a BP 
wave packet and leave a coherent shift in the position of it. This suggests a simple scheme for 
constructing a BP charge qubit to realize a quantum switch, which is capable of controlling the 
coherent transport of one and only one photon in a one-dimensional waveguide. 

PACS numbers: 03.65.Ge, 05.30.Jp, 03.65.Nk, 03.67.a, 



Introduction. Most recently, many theoretical and 
experimental investigations about bound pair (BP) in 
strongly correlated boson systems are carried out [l], 0, 
B HL H 0, HI since the experimental observation of 
atomic BP in optical lattice [1]. Counter intuitively, it 
is found that the trapped rubidium atoms in a three- 
dimensional optical lattice can form a stable BP, even 
though in free space the two atoms would have repelled 
each other. For the problems BP, we can cast back for 
much earlier investigations of ^-pairing states in Hubbard 
model for electrons, which possess off-diagonal long-range 
order (ODLRO) 9]. Actually, the basic physics of both 
fermion and boson BPs is that, the periodic potential 
suppresses the single particle tunnelling across the bar- 
rier, a process that would lead to a decay of the pair. In- 
teresting questions are whether such a BP as composite 
particle will occur in moderate U system and whether 
it can exist stably as a wave packet. It is crucial for 
quantum information processing since the Bose Hubbard 
model is the simplest model capturing the main physics 
of not only cold atoms i n op tical lattice also photons in 
nonlinear waveg uide [id, IHtEl- 



In this paper we will present some exact results ob- 
tained by Bethe ansatz concerning the two-particle prob- 
lem. We demonstrate the existence of a stable BP wave 
packet in Bose Hubbard model with arbitrary on-site in- 
teraction U. It is found that the most statble BP wave 
packets refer to different regions of a bound pair band 
(BPB) and have different group velocities as U varies 
from zero to infinity, but spread to the same fidelity when 
they travel over the same distance. This feature allows 
the BP wave packet as a new object to be a flying and 
stationary qubit in quantum device. We also investigate 
the scattering between a BP wave packet and a single 
particle (SP) in large-C/ limit. It is found that an in- 
coming SP wave packet can coherently pass through a 
BP and leave a coherent shift in the position of the BP, 
which arises from the exotic effective exchange interac- 
tion between them. Furthermore, utilizing on-site U one 
can confine a BP, rather than a SP, in two sites to form 
a charge qubit. This suggests a simple scheme to re- 
alize a quantum switch, which is capable of controlling 



the coherent transport of one and only one photon in a 
one-dimensional waveguide. 

Wave packet in bound-pair band. The simplest model 
capturing some physics of the nonlinearity of photons in 
a coupled cavity array and cold atoms in optical lattice 
is a Bose Hubbard model. The Hamiltonian H is written 
as follows: 
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where a\ is the creation operator of the boson at zth 
site, the tunnelling strength and on-site interaction be- 
tween bosons are denoted by k and U. For the sake of 
clarity and simplicity, we only consider odd-site system 
with N = 2No + 1, and periodic boundary conditions? 
a-N+i = 0,1. 

Consider the two-particle problem, a state in the two- 
particle Hilbert space can be written as 
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where k = 2irn/N, n £ [1,-N] denotes the momentum, 
and re [I, No — 1] is the distance between two particles. 
Due to the translational symmetry of the present system, 
the Schrodinger equations for / fc (r), r G [0,iVo — 1] is 
easily shown to be 

No -i 

[ £ {S^+x + 5 i , r _ 1 )-^ r ,o+T r fc 5 r , J v -£ fc ]/ fe (r) = 0, 

3=0 

(3) 

where = -2v / 2k cos(fc/2) for r = 0, and -2kcos(A;/2) 
for r ^ 0, respectively. Obviously, for an arbitrary fc, the 
solution of (|3|) is equivalent to that of a non-interacting 
iVo-site tight-binding chain with nearest neighbor hop- 
ping amplitude Tj, on-site potentials U and — 2kcos(/c/2) 
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FIG. 1: (Color online) (a) Plots of (rii(i)) for the wave pack- 
ets with a = 2/15 at time t = (dot) and t = t = 15/k 
(empty circle) in the systems with U = 0.1~20 in the unit of 
n. (b) Plots of v g in the unit of k (dashed line) and size Afc 
(solid line) for the BP Gaussian wave packets, (c) Schematic 
illustration of the exchange interaction between SP and BP. 



at two ends respectively. In this work, we focus our study 
on the bound states. In each fc-invariant subspace, there 
exists only one bound state for nonzero U, which can 
be obtained via Bethe ansatz method. And all the N 
bound states, indexed by k, constitute a bound-pair band 
(BPB). 

For a large N system, the pair-bound band can be 
expressed as 



£fc = sgn(J7)\/ U 2 + 16k 2 cos 2 — 



(4) 



with the wave function 
/ fc (r)c[sgn(C fc )r(l + C 2 fe 



1, (r = 0) 

(5) 

where ( k — incos(k/2)/U and \i k — ln[l/£ fe + 



1 + (1/Cfe) ]• The spectrum of BP ((4]) is in agreement 
with that obtained from the Green's function method [l[ . 
The size of the BP for every bound state can be charac- 
terized by 
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which depends not only on k/U but also on fc. It can 
be seen that even for weak U, the size of BP still re- 
mains small for long- wave eigenstates, which is crucial 
for the following discussion. On the other hand, for each 
eigenstate, the BP is delocalized as a composite parti- 
cle. Nevertheless, it has been argued that approximately 



nonspreading wave packet can be achieved by a super- 
position of eigenstates within a linear region, so as to 
the populated energy levels are equally spaced 3, 3j- 
Note that there exists a linear region in the vicinity of 
fco in the BPB spectrum (fj|]) for any value of U. Here ko 
is determined by the condition {d 2 e k j '<9fc 2 ) k _ k = or 

its more explicit form cos fco — V 7 ? 2 — 1 — Vi where rj = 
(U 2 /8k 2 + 1). Within such a region, a Gaussian wave 
packet can be constructed in the form 
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where N c £ [1,A] is the center of it in real space, 



and f2 = J2h 



- i (fc-fe ) 2 



is the normalization factor. 



The dynamics of such a wave packet is governed by the 
effective Hamiltonian H e ff = J2 k ik IV'fc) (V'fcl approxi- 
mately. Here the effective linear dispersion relation is 
Ik = £ft + v g (fc - fco), where 
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is the group velocity of the wave packet ((7j) in real space. 

Now we firstly investigate the dynamics of such wave 
packet for any value of U under the linear approxima- 
tion. Taking the state ([7]) as an initial state = 0) 

its time evolution driven by H e ff presents ^(t) 

e^ H -ff t \<i>(k ,N c )) = e lv \<f>(k ,N c + v g t)). The over- 
all phase factor e lip has no effect on the final result. It is 
obvious that the wave packet moves along the ring with 
velocity v g . In this sense, the time evolution of some 
states governed by H e ff can be described as a spatial 
translation by the operator hi(t) = exr>(—ikv g t) = T(l) 
with a displacement I = v g t. This shows that the shape 
of the wave packet in the real space does not change 
approximately during its travel. However, for the ex- 



act time evolution, state ^(t)) 



-iHt 



$(fc ,Ay> is 



slightly different from state S&(t)) due to the nonlinear- 

ity of the dispersion ((4]) . The overlap between two states 
that evolve from the same initial wave function under 
two different Hamiltonians H and H e ff, respectively, is 
defined as the Loschmidt echo (LE) or quantum fidelity 

F(t) — {^(t) ty(t)\ which can be employed to depict the 
deformation of a travelling wave packet. A straightfor- 
ward calculation shows that 



F(t) 



(fc - fc ) Vgt 



(9) 



which is based on the fact (c> 3 £fc/<9fc 3 ) fc _ fe = — v g . 

Remarkably, the fact that the fidelity (J9j) only depends 
on Vgt means that the wave packets with fixed a but dif- 
ferent fco share the same fidelities after they travel the 
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same distance I = v g t. It indicates that a slower wave 
packet in strong U system has longer life time comparing 
to a faster one in a weak U system. This feature can be 
utilized to quantum information processing: weak U sys- 
tem can be a quantum channel for quantum state trans- 
fer, while strong U system can be employed for quantum 
state storage. 

The above discussion tells us that a state of type ([7|) 
is non-spreading only within certain approximate limits. 
Next we investigate the profile of such a state in real 
space. It is well known that for a SP case, if we replace 

\ip k ) as \k) = l/VNJ2j elkja ) \ Vac )' the SP wave Ac- 
tion is 

|<Mfc , N c )) = 4=£ e -Ttf-".) a +*oi t \ Vac) t (10) 

3 

which is also a wave packet of Gaussian type For 
two-particle state ([7]), its profile in real space can be de- 
scribed by the distribution of the average particle density 

K(t)) = (*(t)|afa i |*(t)) = |a i |*(t))| 2 . (11) 

In Fig. 1 (a) we plot (ni(t)) for the wave packets with 
a = 2/15 at time t = and r = 15/ k in the systems 
with U = 0.1, 0.5, 1, 5, 10, and 20 in the unit of k. 
It shows that the profile of the wave packets in the real 
space is also Gaussian type and non-spreading. It also 
indicates that the shape of the wave packets does not 
change apparently for different U, which is in agreement 
with the following observation from Eq. (fT2|) that the 
size of a BP remains small for \U\ £3 0.2k. 

For a given a, the size of a composite particle in the 
form of a nonspreading wave packet is a function only of 
the ratio k/U 

1 / / 16k 2 , . 

which determines the size of the wave packet. To demon- 
strate the features of a BP wave packet for arbitrary U 
system, its velocity v g and size Afe are plotted in Fig. 
1 (b). For U/k = 0.2 - 10, we have A feo =0.2-3.1, 
which is sufficiently small that we can have many pairs 
in the lattice without having substantial overlap between 
them. It implies a new phase, a gas of BPs which has 
been predicted in large U limits 3], can also exist in 
moderate U system. It is worthy to stress that, although 
a SP wave packet (flQ|) and a BP wave packet ((7]) share 
some common properties, on-site interaction U is able 
to govern a BP rather than a SP wave packet. In this 
sense, the SP and BP wave packets can be regarded 
as two different types of particles. Nevertheless the in- 
teraction between them is exotic since the constituent of 
BP is essentially SP. 

Coherent shift. In the following we restrict ourselves 
to large U limit. Consider a three-body problem. The 
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FIG. 2: (Color online) The stroboscopic picture of the profiles 
of evolving SP and BP wave packets obtained by numerical 
simulations: line a denotes the center of SP wave packet, 
while lines b and c denote the centers of BP wave packets 
before and after scattering. 

spectrum consists of three bands around 0, U and 2U. 
We are interested in the middle band, which corresponds 
to a SP and a BP. Using perturbation method, the cor- 
responding effective Hamiltonian is 

N 2 N 

i=l i=l 

N N 

-V2k &t +1 6j5|Sj+i + h.c. + U*^2 

i=l i=l 

where hi and hi denote the hardcore bosons satisfying 
the following commutation relations [a^* , at] = = 
[bj,b\] = 0, (i ? j); {a h al} = = {Kb}} = 0. 

The first two terms describe the hopping of SP and BP, 
the third term describe the interaction between the two 
kinds of particles, the process of which is schematically 
illustrated in Fig. 1 (c). Now we focus on the scatter- 
ing between SP and BP wave packets. In short time 
duration, a BP is relative stationary comparing with a 
moving SP wave packet. Then the swapping operation 
&J+i&i5|ai+i allows the incoming SP wave packet "pass 
through" the BP and shift its position with a unit lattice 
spacing. To demonstrate this process, numerical simula- 
tion is performed for the time evolution of such two wave 
packets. The initial state is \4>(ir /2, N sp )) |$(tt/2, N bp )} 
with N sp <C Nbp. In the simplest case of replac- 
ing the swapping term by k, the final state is 
\4>(TT/2,N sp ))mir/2,N bp -l)) with N sp > N bp . Actu- 
ally, factor \[2 can cause a slight reflection of the incom- 
ing SP wave packet from the above fact. Fig. 2 is the 
stroboscopic picture of the profiles of two evolving wave 
packets obtained by numerical simulations for the Hamil- 
tonian (|13[) : line a denotes the center of SP wave packet, 
while lines b and c denote the centers of BP wave packets 



before and after scattering. It is clear that the incoming 
SP wave packet keeps the same speed during the whole 
process, while the BP wave packet get a coherent shift 
with a unit of lattice spacing. 

Bound-pair charge qubit. Now we apply the novel fea- 
ture of coherent shift of a BP to construct a quantum 
switch. Considering a 4-site chain in strong on-site inter- 
action limit, the dynamics of a SP and a BP obeys the 
Hamiltonian 



H C q = - 



k (ala s+ i + al +2 a s+3 \ - K al +1 a s+2 (14) 
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which is schematically illustrated in Fig. 3 (a). We fo- 
cus on the case of that there is a single BP in the site 
s + 1 and s + 2, which can be realized under the condi- 
tion U ^> Uq. Notice that this setup is equivalent to the 
system of confining a composite in an effective double- 
well potential and can be regarded as a charge qubit. 
Such a qubit has a novel feature due to the coherent 
shift induced by the scattering with a SP. To demon- 
strate this, we take v2ko = K f° r simplicity and study 
the dynamical process via time evolution. We embed 
such a charge qubit into a chain as illustrated in Fig. 3 
(b, c). Let us firstly assume that initially the qubit is 
in the "right" state \R) = b\ +2 |Vac), while a SP wave 
packet of type dTUJ) |0(vr/2,iV c -< sj) = \<p(n/2,L)) (sim- 
ilarly, we define \4>(±n/2, N c >- s + 2)) = \4>{±ir/2,R))) 
is coming from the left. Comparing to the speed of the 
SP wave packet v g = 2k, state \R) can be regarded as 
a stationary state during the whole scattering process. 
Then according to the Hamilltonian Hcq, the incom- 
ing wave will pass through the qubit freely but leaves 
the qubit to be in the "left" state \L) = b\ +1 \ Vac), 
i.e., \(j)(n/2,L))\R) — ► \4>(n/2,L))\L). In contrast, 
if the qubit is in state \L), the scattering process is 
|</>(7r/2, £)) \L) — ► \4>{--k/2,L))\L), i.e., the incoming 
wave packet is totally reflected and the qubit remains 
to be in state \L). These two processes are illustrated 
schematically in Fig. 3 (b, c). Remarkably, if a sequent 
wave packets scatter with the BP qubit in \R) state, the 
first one can pass freely, but the subsequent ones will be 
reflected totally. This suggests a simple scheme to real- 
ize a quantum switch to control the coherent transport 
of a photon in a one-dimensional waveguide. The pho- 
ton blockade [l5| can be utilized to construct a photon- 
pair qubit in coupled-cavity array. Nevertheless, different 
from schemes in Refs. 16|, [17J , our scheme allows one and 
only one photon passing over the switch. 

Conclusion. In conclusion, we have studied the exis- 
tence of localized BP in Bose Hubbard model with arbi- 
trary on-site interaction U . We have shown that BP wave 
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FIG. 3: (Color online) (a) A confined BP as a charge qubit 
with states \R) and \L). (b, c) A BP qubit as a quantum 
switch to control the transport of a SP wave packet. In the 
case with the qubit in state \R), a moving SP wave packet will 
pass through the qubit freely but leaves the qubit to be in the 
\L) state. In contrast, if the qubit is in state \L), the coming 
wave packet will be totally reflected and remains the qubit to 
be in state \L). 



packets refer to different regimes of a bound pair band 
(BPB) and have different group velocities as U varies 
from zero to infinity, but spread to the same fidelity when 
they travel over the same distance. It proposed a new ob- 
ject to be a flying or stationary qubit in quantum device. 
Furthermore, the coherent shift in large- U system sug- 
gests a BP qubit as a quantum switch embeded in a one- 
dimensional waveguide. Our analysis can be extended to 
a fermion Hubbard system with a minor correction. 
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